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Abstract 

Absorption profile of a four-level ladder atomic system interacting with 
three driving fields is studied perturbatively and analytical results are 
presented. Numerical results where the driving field strengths are treated 
upto all orders are presented. The absorption features is studied in two 
regimes, i) the weak middle transition coupling, i.e. Q2 « ^1,3 and ii) 
the strong middle transition coupling Q% >> ^1,3- In case i), it is shown 
that the ground state absorption and the saturation characteristics of the 
population of level 2 reveal deviation due to the presence of upper level 
couplings. In particular, the saturation curve for the population of level 
2 shows a dip for fii = ^3. While the populations of levels 3 and 4 show 
a maxima when this resonance condition is satisfied. Thus the resonance 
condition provides a criterion for maximally populating the upper levels. 
A second order perturbation calculation reveals the nature of this minima 
(maxima). In the second case, I report two important features: a) Filter- 
ing of the Aulter-Townes doublet in the three-peak absorption profile of 
the ground state, which is achieved by detuning only the upper most cou- 
pling field, and b) control of line- width by controlling the strength of the 
upper coupling fields. This filtering technique coupled with the control of 
linewidth could prove to be very useful for high resolution studies. 



1 Introduction 

It is well known that the control of absorption and emission properties of atoms 
may be achieved by controlling the nature of atomic coherence. Driving field 
strengths and detunings provide the necessary tuning parameters to achieve this 
control. Multi-level atoms have the advantage of providing a larger set of tuning 
parameters. Three level atoms show a class of phenomena which are absorption 
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inhibitive £Q while four-level systems have shown absorption inducive features 
at three photon resonance [2]. The narrow absorption features which have been 
seen in four level systems is attributed to the interaction of double dark res- 
onances 13- Scully et al 0] report the occurrence of Doppler-free absorption, 
which has also been reported earlier |2J. Other features which have been re- 
ported in four-level systems include coherence switching py, photon switching 
(two-photon absorption) Jj], fast switching of nonlinear absorption [7j, enhance- 
ment and suppression of two photon absorption [Hj and the occurrence of three 
peaked absorption [2J El E3 [U\ ■ I n this paper we investigate the absorption 
properties and the modification of the ground state absorption due to the ad- 
ditional upper transition couplings. In section II the details of the model is 
described and the equation of motion is set up. Previous studies have 
revealed that the dynamics of four-level atoms interacting with three driving 
fields may be classified into two broad domains, i) the weak middle transition 
coupling and ii) the strong middle transition coupling. Hence, the analysis is 
presented separately for these two cases in sections III and IV respectively. In 
both the cases analytical results are derived perturbatively and the qualitative 
features are compared with exact numerical calculations which are illustrated 
graphically. 

2 Equation of motion. 

The model we consider consists of a four-level ladder system shown in Fig 1. 
The only dipole allowed transitions are 1 «-» 2, 2 <-> 3, and 3^4 coupled 
respectively by the Rabi frequencies Qi,D,2,^3 and the Bohr frequencies are 
denoted by u>i,i — 1,3. The decay constants corresponding to the respective 
levels are denoted by 1^, i=l,4. The frequency of the applied fields are denoted 
by u>Lt, i = 1) 3 and the detunings are denoted by A<, i — 1,3 corresponding to 
the three couplings. The analysis presented here can also in principle be valid 
to any other type of 4-level system with transitions such that not more than 
two dipole transitions share a single level. This could include 'N' type systems 
as well as the mirror reflected 'N' type systems. This particular ladder system 
may be identified, for example with the Rb hyper fine levels. For instance, 
the 5si,5p3 and 5ds could correspond to the three levels with an additional 

2 2 2 

hyperfine level corresponding to either 5pa or 5ds . In fact the decay parameters 
of the model have been chosen to be the same as the decay constants of these 
levels, T2 = 67, T3 = 7 and T4 = 7 in terms of the ground state life time 7 
which is chosen to be one for convenience. The density matrix equations (|12j) 
are written in the rotating wave approximation as 



dpi2 
dt 



(-iAi - r 2 /2)pi2 - i^l(p22 ^ /Oil) + «^2Pl3 



(-iA 2 - (r 2 + T 3 )/2)p 2S - iflipis - i^2(P33 - P22) + «^3P24 
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Figure 1: Four- level system interacting with three driving fields of strengths 
Oi, f2 2 , ^3 • r 2 , T 3 , T 4 are the decay constants of the corresponding levels. 



dp; 



'34 



dt 

dpi3 
dt 

dt 

dp 2 4 

dt 

dp 2 2 

dt 

dp33 

dt 

dp44 

dt 



= (-iA 3 - (r 3 + r4)/2)/934 - ifi 2 p 24 - i^z{p4i - P33) 

= (-i(Ai + A 2 ) - r 3 /2)pi 3 - iflip 2 3 + i^2Pl2 + «^3Pl4) 

= (-*(Ai + A 2 + A 3 ) - T 4 /2) Pl4 - m lP 24 + i^3p\3 

= {-i(A 2 + A 3 ) - (r 2 + r 4 )/2) ( o 24 - iSlipu - «0 2 p 3 4 + i«3P2 3 ) 

= -r 2 p 22 + ifil(p21 - P12) + «^2(P23 - P32) 

= -r 3 p 33 + if2 3 (p 3 4 - P43) - «^2(P23 - P32) 

= -r 3 p 44 - zfi 3 (p 34 - p 43 ) (1) 



where are the decay constants of the levels i and Aj = u)%,i+\ — Wi are the 
laser detunings (pij = pj* and Trp =1). We study the steady state behaviour 
of the system for various values of the driving field strengths and detunings. 
In the general situation, we present the numerical solutions to the populations 
and coherences where the driving field strengths are treated up to all orders. 
The modification of the absorption profile of the ground state as well as the 
upper levels will be studied analytically by perturbatively solving for the density 
matrix elements in Q 2 or fii as the case may be. 
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3 Weak middle transition coupling. 



In the case when the transition 2 <-> 3 is weakly coupled, i.e. when f2 2 is 
very small compared to f2 1; f2 3 , the equation of motion for the density matrix 
elements can be solved perturbatively in f2 2 ■ The set of fifteen coupled equations 
decouples into three subsystems of equations which are solved independently 
to obtain analytical solutions. In fact, the density matrix elements given by 
the sets {pn, p 22 , P12, P2i},{p33, P34, P43, P44} and {p 23 , Pi3, Pu, P24} form three 
subsystems. The zeroth and first order solutions for pij,i,j = 3,4 remain zero, 
while the coherences {P23, P24, P13, P14} are nonzero in the first order and are 
listed in the Appendix. The steady state second order solutions for py are listed 
below: 

P$ = (2f 2 n?(i - P $ P i 2 4 } ) - n 2 Imp^ (A 2 + f D + 

2ft 1 r* 2 (A 1 7TOp( 1 3 ) + f 2 Rep[\ ) )/V 1 
pf 2 } = -i(f 2 ni(iAi - f 2 )(l - p$ - pf}) + f 2 fi 2 (iAi - f 2 )pW 

+2fiifi 2 (*Ai - f 2 )/mp 2 1 3 ) - 2f7 2 fi 2 (p^ 3 ) - p i ^)/V 1 
pf 3 } = (r 4 2 (A 2 + (f 3 + f 4 ) 2 ) - 2fi 2 r» 2 (f 3 + f 4 ))/mp 2 1 3 ) 

+2f! 2 fl 3 r 4 (A 3 /mp 24 ) + (f 3 + T 4 )Rep$)/V 3 
p^l = (2Q 2 n 3 f 3 (A 3 /mp 2 1 4 ) + (f 3 + f 4 )i?ep 2 1 4 ) ) + f> 2 fi 2 (r 3 + f 4 )/mp 2 1 3 ) )/I? 3 
Pa? = 2(2zfi 2 r! 2 (f 4 + f 3 )/mp 24 ) - f 3 f 4 r» 2 (zA 3 - f 3 - f 4 )p 2 1 4 ) 

+f 4 r» 2 r! 3 (zA 3 - (f 3 + f 4 ))/mp 2 1 3 ) )/I? 3 (2) 

where V 1 = f 2 (A 2 + f 2 + 4ft 2 ) and V 3 = 2ft 2 (f 3 + f 4 ) 2 - f 3 f 4 (A 2 + (f 3 + 
r 4 ) 2 ). Here, for convenience we have relabeled Tj/2 everywhere by I\ . The 
second order steady state solutions for the set of coherences {P23, P13, P24, P14} 
are the same as the first order solutions given in the Appendix. 

Note that, in the limit £! 2 — > and ft 3 — > 0, the solution for pi 2 and 
P22 are given by the zeroth order solutions which is the usual absorption and 
population of two-level system driven by a field of strength Oi. There is very 
little modification in pi 2 and p 22 due to the first order contribution. However the 
inclusion of the second order contribution introduces additional terms. Before 
investigating this further, let us look at the graphical illustration of the ground 
state absorption as a function of detuning which has been solved numerically for 
all orders in the driving field strengths. In Fig 2 the curve corresponding to (a) 
shows the absorption in the presence of only 1 ^> 2 coupling, (b) corresponds 
to the absorption of the ground state (g.s.) in the presence of the 1 <-+ 2 as well 
as the 2^3 couplings and (c) shows the modification of the g.s. absorption in 
the presence of the 1 <-> 2 the first excited state coupling 2^3 and the 3^4 
coupling. Here the small dip in the line center is not due to population trapping 
but due to the population transfer to the upper levels which is because of the 
inclusion of the upper most transition coupling. This becomes more clear by 
looking at equations (2) and (3) which will be discussed shortly. Let us look at 
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the variation of the populations of various levels as a function of fii and for 
very weak f^- Fig3-5 show the variation of populations w.r.t. Oi and O3 which 
have been obtained numerically. Note that P22 shows a dip in the absorption 
when f2i = Q3. There is however a corresponding increase in the populations of 
the upper levels. Fig4 and Fig5 show a sharp increase in the absorption by level 
3 and 4 respectively as fii approaches £! 3 and shows a maxima for Six = fl 3 
when f2 2 — 2. 

One can explain this analytically by looking at the second order solutions. 
Consider the situation when all the detunings are zero and ignoring the contri- 
bution of the decay parameters 0(Tf) i.e. assuming Q ly3 » r 2 ,3,4 for A; = 0, 
the populations reduce to the simple form 



(2) _ (0) nfofolcFsp^-niJmp^) 

^22 ~ P22 ' f 4 (o2-o2)+g 

Tm J 2 ) ~ TmnW , 4O^2l» 1 ((r 2 p< 2 ) (2r3-r4)+/mp<" 1 ) (-2r 2 +r 4 )o 1 

lmp 21 ~ lmp 21 H r 4 (a?-a|)+6 

(2) _ r 4 o^o 3 (-P 2 2 ) r 3 +niJmpW) 
irnp M ~ 2r 4 n§(n?-n§)+c; 

(2) ^ f2 2 (2H| (F 4 -r 3 )p^° ) +2n§ fmp|j ) 

^33 ~ 2r 4 n^(nj-n^)+e 

(2) _ 2Q^^(r 3 pW-0 1 JmpW) 
^44 ~ 2r 4 02(02_ 2 )+e 



(3) 



where = 2fi i/( r 2 +_4«?) and /m^"' = r 2 Qi/(r| + n\) for Ai = and 
Q includes terms of the 0(T 2 ). It is obvious from these expressions that there 
is a resonance at Hi = fl^. The populations p^ and p$ show a sharp rise as 
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Figure 4: p S3 as a function of the driving field strengths f^i and f2 3 with £1 2 = 2. 
Occurence of a maxima for Oi = f^. 
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Figure 5: P44 as a function of the driving field strengths f2 1 and f2 3 with fi 2 
and T 2 = 6,r 3 = T 4 = 1. Occurence of a maxima for Sli = Q 3 . 
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Q. 

Figure 7: P33 as a function of the driving field strengths fl a) Vl 2 = 2, Hi = 

20,o = n 3 b) n 3 = 20, n 2 = 2,n = Oi c) o 2 = 8,Oi = 20,0 = n 3 d) 

3 = 20, n 2 = 8, n = fii. Here, T 2 = 6, T 3 = T 4 = 1. 




Figure 8: P44 as a function of the driving field strengths Oi. a) Sl 2 = 2, Oi = 
20, Q = Q 3 b) Q 3 = 20,^2 = 2,fi = Oi c) 2 = 8,fii = 20, Q = Q 3 d) 
3 = 20, n 2 = 8, n = S7i. Here, T 2 = 6, T 3 = T 4 = 1. 
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fil — > since the numerator is positive. Thus to populate the levels 3 and 4 

maximally when the middle level coupling is weak, it is sufficient to match the 

(2) 

requirement f^i = f2 3 . The absorption profile imp 34 also shows a similar trend. 

(2) (2) 

While the second term in p 22 and p 21 are negative and the role of the resonance 
in the denominator is to introduce a minima when Qi = Q3. Thus the role of 
the second order term in this case, is to introduce a dip in the line center in an 
otherwise two-level absorption profile and hence a depletion in the population 
of level 2. As one increases J7 2 this perturbative expansion is no longer valid. 
Moreover, the numerical calculations upto all orders in fi 2 reveal the fact that 
the resonance condition becomes weaker and the sharp features vanish. This is 
demonstrated graphically in Fig6-8. 

Thus, in a four- level system with middle transition coupling very weak, the 
condition for maximally populating the upper levels is to satisfy this reso- 
nance condition. The population transfer to higher levels occurs at the cost 
of depletion of population in the second level. This modification is mediated 
by the contribution of the coherences, which at three photon resonance are 
given by p^ « fl 2 Imp^{ /Q , P13 ~ n^p^ /Go, P2I ~ ^2^3/4 2 } /So where 

So = r 4 (fi? - ft§). 



4 Strong middle transition coupling 

I next consider the case when the 2 — > 3 transition is strongly coupled i.e., 
Q2 » f2i,f23. Treating Qi perturbatively, the first order absorption of the 
ground state takes the simple form 

(i) _ L 12 + »|/L 123 0i 
lm P2i ~ (^ + L 1 £ 12 + n§L 1 /i 12a ) 1 j 

where L x = (-iAi - T 2 /2),L 12 = H(A 2 + Ai) - T 3 /2) and L 123 = (-i(Ai + 
A 2 + A3) — T4/2). We now look at some simple situations, for example when 
£1 2 —* 0, fls — > 0, the absorption becomes proportional to r 2 J7i/(Af + r 2 ) which 
is the ground state absorption for a weak excitation and observe that the width 
is proportional to T 2 . Next, consider only — * 0. The absorption is now given 

by 

T = »i(Afr 2 + r 3 (r 3 r 2 + »l)) 
921 Af + Af(fl + r|-2n2) + (r 2 r 3 + n|) 2 { > 

where A 2 is assumed to be zero. The denominator is quadratic in A^ and hence 
has two peaks corresponding to the usual Autler-Townes doublet. The width 
of this EIT window is proportional to ^/(f \ + f § - 2S^|) 2 - 4(f 2 f 3 + Q|). The 
absorption for Ai = (at the line center) is proportional to r2if3/(f 2 f3 + f2|) 
which is negligible since it is of the order Vli/Vt^. The absorption in the presence 
of the third driving field for A 2 = 0, A3 = takes the simplified form 

j (i) _ fli(Aj(Ajf 2 + f 2 (f I + f I) + f 3 nj - gTgOf) + Ti) 
mp21 A?(Af - T p - Q? 2 - tti) 2 + (A\T S - f 2 f 3 f 4 - f 4 n| - f 2 r>|) 2 1 ' 
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where T s = f 2 + f 3 + f 4 , T p = f 2 f 3 + f 3 f 4 + f 4 f 2 and Ti = f 2 f§f| + 
r 3 rffi| + 2r 2 r 3 r 4 fi§ + r 4 Sl§il| + r 2 fi^). The denominator is a cubic equation 
in Af with the three roots corresponding to the three absorption peaks, two 
of them correspond to the usual Autler -Townes doublet and the absorption 
at the line center is due to the three photon interaction and is proportional to 
f2if2§/ (f^fi! + ^2^3) which occurs when all the detunings are zero. This three 
peaked absorption in four level systems has been observed earlier and has been 
frequently referred to as the splitting of the EIT window. While the role of the 
strong 2 «-> 3 coupling was to introduce transparency at the line center, the 
role of the 3 <-> 4 coupling is to induce a narrow absorption at the line center 
within the EIT window. Thus the two photon absorption was responsible for 
the coherence p 34 which in turn lead to non-absorption while, the three photon 
absorption which introduced a coherence between levels 1 and 4 (pn) which in 
turn induced this sharp absorption feature. Thus there is a contrast in the nature 
of the quantum interference in the case of two and three photon interaction. 

On the other hand consider the case when the detuning A 3 is nonzero. Let 
Ai = 0, A 2 = 0, the absorption now reduces to 



(1 ) _ Oi(A|f 3 (f 2 f 3 + flf) + (f 3 f 4 + n|)(f 4 »1 + f 2 (f 3 f 4 + »§))) , 

mp21 Ai(f 2 f 3 + ni)2 + (r 4 f)i + f 2 (f 3 f4 + fi|)) 2 [) 

This is a Lorentzian with width proportional to F 2 f2 3 /f2 2 . This implies that the 
line width can be subnatural (< T 2 ) for f2 3 < fi 2 . Notice the novel feature here 
which is the absence of the Autler- Townes doublet unlike the case when only 
Ai was nonzero. Thus, we now have obtained a powerful technique for filtering 
the Autler- Townes doublet and retaining only the narrow absorption at the line 
center in addition to gaining control over the linewidth. This could prove to be 
a useful technique in high resolution spectroscopy. 

A numerical illustration of these novel features is presented in Fig 9 where the 
steady state solutions are obtained for all orders in the driving field strengths. 
Figure 9 clearly shows the difference in the behaviour depending on the strength 
of the middle transition coupling, a) corresponds to weak f2 2 while b) corre- 
sponds to strong 0,2 . Fig 9 A-b shows the three peaked absorption while Fig9B-b 
shows the elimination of the Autler- Townes doublet by detuning the driving field 
coupling 3-4 transition as predicted by the perturbative calculations. Likewise 
for small f2 2 , only the narrow absorption dip at the line center is seen by de- 
tuning A 3 (Fig 9A-a and Fig9B-a). Again, Imp 32 for small 17 2 is very strong 
(Fig9C-a) as the perturbation calculation indicate (section III) and this absorp- 
tion is not very significant for large J7 2 since the resonance condition weakens 
and the contribution of higher order terms in f2 2 dominate. Lastly, the absorp- 
tion profile of the uppermost transition does not seem to vary much with cither 
of the detunings Ai or A 3 (Fig9-Eb,Fig9F-b). This is very much similar to 
the narrow g.s. absorption at the line center. Further, both the absorption of 
the g.s. (line center) and the upper most transition 7m/? 43 are unaffected by the 
detunings A1/A3 since they occur at exact three photon resonance. One would 
expect this absorption feature also to be Doppler-free like the g.s. absorption 

El 
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Figure 9: The variation of absorption for (a)02 = 2, Oi = O3 = 20 and (b) 
0,2 = 20, 51i = O3 = 4 .A)Imp2i with Ai B)7mp2i with ^.Gjlmp^ with Ai 
D)/mp32 with A 3 E)7mp43 with Ai F)/mp4 3 with A 3 with the Ti,i = 2,4 the 
same as in the previous figures. 
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Figure 12: Imp4 3 as a function of ^3 and A 3 for parameter values the same as 
in previous figure. 

Next the variation of the absorption of the g.s., Imp 2 i, the first excited state 
Imp32 and the second excited state Imp^ with the driving field strength f2 3 
and the detuning A3 are numerically illustrated in FiglO-12. The purpose here, 
is to demonstrate the three photon effect. While the ground state absorption 
increases with fi 3 , the middle and the upper transition absorptions rise sharply 
for small Q3 and fall rapidly to zero. All the absorption occurs only for A3 = 0. 
The variation with A 2 shows lot more structure and will be more appropriate 
to present it else where. 

To summarise, it is shown that by satisfying the resonance condition Oi = 
£! 3 one can maximally populate the upper levels at the cost of depleting the 
population of level 2 (for small f^)- In the case of large Sl 2 it is shown that one 
can get narrow absorption features at the line center by filtering of the Autler- 
Townes doublet. One also has a control over the line- width by controlling Q 2 
and Q3. The qualitative features predicted by perturbative calculations seem 
to agree very well with the numerical results. 

5 Appendix 

The zeroth order equation in Q 2 are given by 



dt 
dt 



12 



- (-*Ai - r 2 /2)p 



(0) 



'12 
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g t = A 1 - r 2/2)p21 + ^l(P 2 2 - Pll ) 

the rest of the density matrix elements are zero (steady state). The first order 
equation do not contribute to the solutions of Pij,i,j = 1,2. The density matrix 
elements Pij,i,j = 3,4 remain zero. Hence we furnish below the first order 
equations for the coherences only. 

8 (1) 

= HA 2 -(r 2 + r 3 )/2)p 23 -;Q lP « -in 2 (- P %) + in 3P $ 

8 (1) 

= H(Ax + A 2 ) - T 3 /2)p^ - ifi lP $ + ifi 2 p<°> + ifisp^) 

Q (1) 

- (-i(Ai + A 2 + A 3 ) - r 4 /2)p« - ifi lA &> + ifiap^ 

Q (1) 

- H(A 2 + A 3 )-(r 2 + r 4 )/2)p 2 1 4 ) - J ^ 1 4 ) + *f7 3 p 2 1 3 ) ) 



(8) 



the first order steady state solutions of these equations are given by 

P g ] = n 2 (Q lP ^(nl -n 2 + l 123 l 23 ) + 

i P ( 22 \l 12 l 23 l 123 + L 12 n 2 + L 23 nl))/v 2 
P $ = n 2 {i P ( £{L 123 L 2 L 23 + L 2 nj + L 123 nj) - 

n lP <£(nl-nl + L 123 L 23 ))/v 2 
til = -^Mp ( 22 ] (l 12 l 123 + n 2 ) + 

in lP ^(L 123 + L 2 ))/V 2 

til = n 2 n 3 (ti 2 \nl -n 2 - l 12 l 123 ) + 

in lP ^(L 123 + L 2 ))/V 2 

where V 2 = ((ft 2 - ft 2 ) 2 + Vl\{L 123 L 23 + L 12 L 2 ) + Q 2 3 (L 12 L 123 + L 2 L 23 ) + 
L 12 L 23 L 123 L 2 ) andL 2 = (-iA 2 -(r 2 +r 3 )/2), L 23 - H(A 2 +A 3 )-(r 2 +r 4 )/2) 
while L12 and Li 23 are given in the text. The second order equations are given 
by 

8 (2) 

-g 2 - = (-iA 1 -T 2 /2)p^-iQ 1 (p^-til) + in 2 ti 3 ) 

r, (2) 

°^22 _ p (2) - / (2) (2)x . ,(1) 

— -i2P 22 +"'HP 2 1 - Pl 2 J + «' 2 (,P 23 - P 32 j 

-gi- = (-zA 3 - (T 3 + r 4 )/2)4 2 4 } - ifi 2 p£ - ^(^ - p 3 2) ) 

a (2) 

= -T^+m^-p^-i^til-ti^) 

a (2) 

^ = -T 3 ^-in 3 (til-p^) 
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So far, we dealt with equations where was treated perturbatively. I furnish 
below the relevant first order equations where f2i is treated perturbatively, 

^ = (-iAi-ra^p^-ifii+inap^ 

a (1) 

^ = (-i(Ax + A 2 ) - r 3 /2)p« + zf7 2 p« + ifigpft) 

X) (1) 

= (-i(A 1 + A 2 + A 3 )-r 4 /2)pft>+tfi3pi 1 3 ) 
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